Universal magnetic properties of frustrated quantum 
antiferromagnets in two dimensions 

Andrey V. Chubukov^'^, T. SenthiP and Subir Sachdev^ 
^Departments of Physics and Applied Physics, P.O. Box 208284, 
Yale University, New Haven, CT 06520-8284 
and ^P.L. Kapitza Institute for Physical Problems, Moscow, Russia 

(February 1, 2008) 



Abstract 

We present a theory of frustrated, two-dimensional, quantum antiferromag- 
nets in the vicinity of a quantum transition from a non-cohinear, magneticahy- 
ordered ground state to a quantum-disordered phase. Using a sigma-model 
for bosonic, spin-1/2, spinon fields, we obtain universal scaling forms for a 
variety of observables. Our results are compared with numerical data on the 
spin-1/2 triangular antiferromagnet. 
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A useful classification of two-dimensional, quantum, Heisenberg antiferromagnets is pro- 
vided by the structure of the magnetically-ordered ground state: the spin-condensates on 
the sites can either be coUinear or non-coUinear to each other. CoUinear magnets have been 
extensively studied in recent years and many of their properties are reasonably well under- 
stood. They possess an 0(3)/0(2) order parameter whose fluctuations describe the low 
temperature (T) properties of the magnetically-ordered state The quantum-disordered 
state has only integer spin excitations (the spinous are confined) and spin-Peierls order is 
expected for certain values of the single-site spin p|. The finite-T crossover between these 
two states has also been studied in some detail 0. 

Less is known, however, about non-collinear antiferromagnets, which are the subject 
of this paper. Examples include the triangular, kagome, and square (with first, second, 
and third neighbor interactions) lattices. The magnetically-ordered state completely breaks 
the spin-rotation symmetry, yielding an S0{3) order parameter [Q. Space and time de- 
pendent twists of this order parameter then define three independent spin-stiffnesses, spin- 
susceptibilities, and associated spin-wave velocities. For simplicity, we will restrict our at- 
tention here to magnets with coplanar spins and an internal symmetry (a C^y symmetry 
on the triangular and kagome lattices, and a screw axis symmetry for the incommensurate 
planar spirals on the square lattice), which leads to just two independent stiffnesses {p±, p\\), 
susceptibilities (x±, X\\)^ and spin-wave velocities {c± = {p±/x±Y^'^i cy = (P||/X||)^^^); more 
complicated non-collinear magnets will have similar properties. The long-wavelength action 
for the SO (3) order parameter has an 5*0(3) x 0(2) symmetry, the 0(2) being a continuum 
manifestation of the internal symmetry noted above Q]. A spacetime dimension D = 2 + e 
study of small fluctuations of the SO (3) order parameter about the magnetically-ordered 
state was performed by Azaria et. al. they found that the stiffnesses and suscepti- 
bilities became asymptotically equal upon approaching the critical point separating the 
magnetically-ordered and quantum-disordered phases, with the critical theory possessing an 
enlarged 0(4) symmetry. A large theory based upon Sp{N) symmetry found a sim- 
ilar magnetically-ordered state, but was also able to access the quantum-disordered phase. 
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The latter state was predicted to be a featureless, fully gapped spin-fluid, with unconfined, 
bosonic spin-1/2 spinon excitations. We also note that there are alternative approaches to 
the quantum disordered phase [0 which are quite disconnected from the structure of the 
ordered state. 

In this paper, we shall present a theory of the universal, finite-T properties of non- 
collinear antiferromagnets in the vicinity of the critical point. We will describe the crossover 
from the magnetically-ordered state, with its low-lying spin-wave excitations, to the fully 
gapped quantum-disordered state via an intermediate quantum-critical region. Our results 
are in complete agreement with some previous studies of the magnetically-ordered state 
and the quantum disordered state 0, and establish a fundamental connection between the 
0(4)-symmetric critical point of Ref. and the deconfined bosonic spinous of Ref. 0; a 
related connection was noted recently in Ref. 0]. We will also obtain new results for the 
low T behavior of the dynamic structure factor and uniform susceptibility of magnetically- 
ordered antiferromagnets. 

Our motivation for this study is similar to that for the anolagous recent study of collinear 
antiferromagnets ^ . A given S = 1/2 antiferromagnet may be either magnetically-ordered 
(as is expected for the triangular lattice) or quantum-disordered (the kagome lattice) 0. 
At low T, the magnetically ordered magnet has thermally-excited classical spin-wave fluctu- 
ations (the renormalized-classical (RC) region), while the quantum-disordered magnet has 
only activated deviations from its ground-state properties. At higher T however both these 
magnets are expected to crossover to a quantum-critical [0] (QC) region where classical and 
thermal fluctuations are equally important. Many properties of this region are universal, and 
are thus amenable to numerical and experimental tests. In particular, there are significant 
quantitative differences between the QC behavior of collinear and non-coUinear magnets, 
which are a direct consequence of the presence of deconfined spinous in the latter. 

We begin by presenting our effective action. We choose to describe the local spin config- 
uration by an SU (2) rotation about a reference ordered state. The choice of SU (2) rather 



than S0{3) is significant, and has the immediate consequence of suppressing the vortices ITO 



associated with 7ri{SO{3)) = Z2 for which the SU{2) field is double-valued. This choice is 
motivated partly by the results of Ref. ||^, where vortices were suppressed in the quantum- 
disordered phase by a Higgs condensate. We parametrize the SU (2) matrix by two complex 
numbers zi, Z2 with \zi\'^ + |z2p = 1, and write down the most general, long-wavelength 
action with an SU{2) x 0(2) invariance: 



S = I (fxdr — 



df^z^d^z - ^ [z'^d^z - d^z^ 



(1) 



It is easy to show that g, = l/2p^_^,gr = l/2x![,7x. = (pj - P±)/p±,lr = (xj - X±)/x±, 
where the superscript denotes bare values; note that if the 7^ = 0, 5 has an enlarged 
0(4) symmetry. The action S can be explicitly derived by a long- wavelength analysis of the 



models of Refs. M and B; we have also learned of a recent study of S by Azaria et. al. [11 



The staggered spin-structure factor (wavevectors measured as deviations from the ordering 
wavevector G) can be shown to be the Fourier transform of Re [z\xi,ti)z(x2,T2)^ . Note 
that this is quartic in the consistent with the identification of the z quanta as spin-1/2 
bosonic spinous. 

We studied S by generalizing z to an A^-component, unit-length, complex vector, and 
performing a 1/A^ expansion; S then has a SU{N) x 0(2) invariance, while for 7^ = it is 
invariant under 0(2A^). This method allows us to work directly in D = 2 + 1 and access 
both the QC and RC regions. Note that the extension to large N is different from that used 
in Refs. ||,|12|. 

We expect that S possesses quantum-disordered and magnetically-ordered (with the z 
quanta condensed) as the couplings (say g^) are varied. A key property of the present 
large A^ expansion is that the long-distance physics at the critical point at gx = gc is 
0(2A^)-symmetric. This is manifested in the magnetically ordered phase [g^ < gc) by the 
critical behavior of the stiffnesses. Josephson scaling is obeyed by the fully renormalized 
P\\, Pl, Xh Xi- all of which vanish as {gc — gx)'^, where u is the correlation length exponent 
(z/ = 1 — IG/SNtt"^ + Oi^l/N"^)). However, the relative differences between the stiffnesses also 
vanish at the critical point: we defined Ai = (py — p_l)/ p±, A2 = {x\\ ~ X±)/XJ-y and found 
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Ai=7i(6)-*^+72(6)-^^ 

A2=7i(0)"*^- 272(0)-'^^ (2) 

where 71 = (27^. + 7r)/3; 72 = (7^ — 7t)/3, and C,j is the Josephson length measured in 
lattice units. The positive crossover exponents 0i,2 measure the irrelevancy of the 7^ terms 
in S; the 7^ are actually 'dangerously '-irrelevant as Ai 2 control long-wavelength physics for 
< Qc. To order 1/iV, we found 0i = 1 + 32/37r2A^, 02 = 1 + lU/lhix'^N 0. 
We now present our scaling results for the wavevector- {k) and frequency- {uj) dependent 
staggered (xs) and uniform spin susceptibilities in the vicinity of = Qc- We restrict 



ourselves to Qx < Qc, although more complete results have been obtained WM. We found 



X«(A;,a;) = I 1^ I /c^T (fc, cJ, x, Ai, A2) (3) 



where Nq is the on-site magnetization at T = 0, $!„ are universal functions of the 
dimensional variables k = hc±k/kBT, to = huj/ksT, x = NkBT/4:7ipj_. We found the 
exponent fj = 1 + 32/3tt'^N. The prefactor of $s remains non-singular at Qx = Qc as A^o ~ 
{9c — QxY with 2(3 = (1 + f/)i/. All scaling functions are defined such that they remain finite 
as X — s> 00. As before P], the argument x determines whether the system is in the QC 
(a; ^ 1) or RC (x ^ 1) region. 

An important difference in the above scaling forms from those for collinear magnets 
is in the value of f/. Here we have f] close to unity, while the analogous exponent for collinear 
magnets was close to zero. This is a consequence of the presence here of deconfined spinous: 
it is the z quanta which behave like almost free particles (at T = 0, (z'^z) ~ l/k"^^^ with 77 
close to 0) while the staggered susceptibility is a correlator of a composite operator of two 
spinous [xs ~ l/fc^"'' with f] close to 1). 

We have computed in a expansion to linear order in Ai_2. We describe our 

results as they relate to various observables. 

Correlation length. As in collinear magnets, we define the correlation length, ^, from 



the long-distance e"''/^ decay of the equal-time spin-spin correlation function. We found 
that, to order there is a simple relationship between the values of ^ for collinear and 
non-collinear magnets. For all values of x, the non-collinear C, is precisely 1/2 the previously 
computed C, [@] for the isotropic 0{2N) sigma model. The factor of 1/2 is a signature of 
deconfined spinous. The collinear expression for ^ however must be used with the effective 
values Ps = P±{1 + NAi/{2N^ - 2)), x = + iVA2/(2A^2 _ 2) and c = (ps/xY^^; notice 
also the factor of 4 difference in the couling constant in (|1|) and in 0]. For the physical 
case = 2, we have to first order in Ai 2 that ps = (2p_|_ + p||)/3, c = (2cj_ + c\\)/3, and our 
result for ^ is then consistent in the RC region with that of Azaria et.al. [Q. 

Static uniform susceptibility. The result for Xu is obtained by evaluating the response to 
a vector potential coupled to the conserved charge of the SU{N) symmetry. 

In the RC region [NksT <^ 47rps) "we obtained 



Xu 



(9j^y ( {Nx± + X\\) 
\ h J [ {N+l)x± 



2x± ksT 



(4) 



' A^ 27rc2 

It is worth emphasizing that although we are considering an essentially classical regime, 
the T dependence of Xu is a purely quantum effect - it disappears if the spin-waves had a 
classical, thermal distribution. 

In the QC region [NksT ^ '^T^Ps)-, we found to order 1/A^ 

o.3r 



Xu 



I its:)' k,T^ 

\ he J 47r 



A^ 



+ a x"^/" + . . . 



(5) 



where 9 = 2 log[(V5 + l)/2], x = NkBT/Aixps and a = 0.8 + 0{l/N). Note that the slope 
of the linear in T term is precisely 1/2 of that in the 0{2N) sigma-model 0. The factor of 
1/2 is again a signature of spin-1/2 spinous and should be amenable to experimental tests. 

Staggered dynamic susceptibility and structure factor. In the RC region, the scaling form 
(0) for Xs collapses into a reduced scaling form in which the physical ^, rather than c/ksT 
is the most important length scale ||1|,0]. We found 

■A:^r(Ar-l)n(^+^)/(^-^) 



Xs{k,iuj„ 



Ps{N - 1) [ 47rp, 
e fik^^to^^c) 



X 



(6) 



were / is a scaling function. Note that computations were in fact done only to order 1/N - the 
form at arbitrary follows from a reasonable guess about the wavefunction renormalization 
of the composite field. The overall factor in is chosen such that /(0,0) = 1 + 0{1/N), 
The behavior of f{x,y) at intermediate x,y = 0{1) is rather comphcated, chiefly because 
spin-wave velocity also acquires a substantial downturn renormalization at kC, = 0{1) [|T|. 
However at kC, ~ cu^/c ^ 1, velocity renormalization is irrelevant and we obtained 

fN - 



log(x"^ + y"^) 



(7) 



\N+lJ x2+?/2 

It is not difficult to demonstrate that this result for /(x, y) yields a Xs{k, i^) which is precisely 
the rotationally-averaged spin- wave result for the ordered antiferromagnet, as it of course 
should be at /c^ ^ 1 but kC,j <^ 1. 

We also computed lmxs{k,u!) for real u. In the RC region, we describe the results using 
the dynamic structure factor S{k, u) which satisfies 

where S is straightforwardly related to $5 introduced in (j^). For experimental compar- 
isons, it is sufficient to consider the frequency range u < c/^. We then found cj) = 
uj/27Tck^ for ck ^ u, (in this region of k, collisionless Landau damping is dominant), and 
'B{k,uj) oc (uj^^/c) {{N — l)kBT / A-iTps)^^~^^^'^^'^~^^ for ck ^ uj (the dominant contribution 
is the damping of quasiparticles) . 

Now the QC region. Here we restrict our results to the critical point x = 00, and 
neghgible anisotropy (Ai 2 = 0). For Tick, Tioj ^ fc^r we obtained 

= — (9) 

16(F- (07 + ^(5)2)1-^/2 

where An = I + 0{l/N). At small k and oj, we have Re$, = {V^ / IQtiQ)[1~ (fc^(l + 
2G/v^)— cJ^)/1202 +...] where . . . stand for higher powers of fc, cJ and for regular correc- 
tions in For Im^^ we obtained the following asymptotic limits for large N 
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Im$, 



A^sin(7rr//2) Oju^ - k') _ 
77 =2 ; > 1 

16_ (p2 _ 

ZJ 



cJ< 1, A; > 1 



where 6'(a;) is a step function. In both cases, Landau damping is dominant. Finally, when 
both a; ^ 1 and A; <^ 1, quasiparticle excitations are overdamped and we only know that 
Im$s oc uj. 

In the T = quantum-disordered phase, Im$s shows a clear signature of deconfined 
spinous - the spectral weight at fixed A; is a broadband continuum rather than the delta- 



function peak present in collinear magnets 



Local susceptibility and spin-lattice relaxation. The local dynamic structure factor S{uj) 
is given by S^u) = J d?k S'(A;, co')/47r^. Simple inspection then shows that for u ~ c^^"*^, the 
integration over momentum is also confined to A; ~ and therefore S{uj) is a universal 
observable. The small frequency limit of S{uj) is directly related to the spin-lattice relaxation 
rate of nuclear spins coupled to electronic spins in the antiferromagnet: 1/Ti oc S{uj 0). 
In the RC region, we find using our previous results for the scaling functions that for uo ~ c^^^ 

« M ((i^^) (10) 

For = 2, we then have 1/Ti oc T^/^ ^. 
Deep in the QC region, we found 

where Kiuj) = ZJ sm(7if]/2) /32n at cJ ^ 1, and K(p) = ZJ C^iV^ — l)/647r at cl; -C 1. 
The factors Bn and Cn both behave as 1 + 0{1/N). Clearly then, 1/Ti oc T^. 

Static structure factor. Unlike collinear magnets, the static structure factor ^(A;) = 
/ duj S{k,uj) /2tt in non-coUinear antiferromagnets is nonuniversal because the frequency 
integral over quantum fluctuations is divergent. This follows from the behavior at large 
frequencies where S{k,uj) oc l/uj"^"^ and f/ > 1. The nonuniversality is however more 
relevant for the QC region, where T is the only scale for fluctuations; in the RC region, ^ is 
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exponentially large, and there is a universal contribution to S{k) from classical fluctuations 
which scales as In the RC region we then have S{k) ~ kBTxs{k,0), where Xs{k,0) 
is given by §). At A; = this yields S{0) oc T^n/{n-i) iV = 2, we then obtain 

.5(0) oc T^e- 

Application to the 5 = 1/2 triangular antiferromagnet. We performed a l/S" expansion 
on this antiferromagnet to obtain the T = values of p_L, py, xx? X\\ (all to order l/S"), and 
Nq (to order l/S"^). For S = 1/2 this gave us No = 0.266; x± = 0.07/ Ja^; x = 0.077/ Ja^; 
Ps = 0.087 J; c = {ps/xY^'^ = 1. 06 Ja. For the uniform susceptibility in the RC regime we 
then obtained Xu = {gpB/haf [0.08/J+ O.OSAjijT/ 7^+ 0{T^ / J^)]. On the other hand, in 
the QC regime, we have Xu = {gps/haf [O.lSfc^T/ 7^+ 0.07/ J {keT /2Tip,f-^/'' + ...]. The 
temperature dependence in the subleading term is likely to be quite small in the region of 
experimental interest {ksT ~ 2ttPs), and we can well approximate this term by a constant. 
Note however, that the factor 0.07 is an = oo result - the 1/A^ corrections to this factor 
have not been computed. Further, the correlation length behaves in the RC regime as 
^ ^ 0.25 {Anps/ksTY^'^ explAnps/ksT] where 47rps 1.09J, and deep in the QC region as 
^ = 0.53Ja/kBT. 

Consider now the numerical results for Xu- The data of recent series expansion studies |T3 



show that Xu obeys a Curie- Weiss law at high T, passes through a maximum at T 0.4J, 
and then falls down. The region below the maximum is quite small; nevertheless, we fltted 
this data by a straight line and found 0.13±0.03 for the slope and about 0.06 for the intercept 
- the results are in better agreement with our QC rather than the RC result. Finally, at 
very low T, we expect a crossover to the RC regime, and the corresponding value of Xu at 
T = is also consistent with the data. We also compared the data for the correlation length 
and 5'(0) at ksT ~ 0.4J and found rough consistency with our expressions in the crossover 
region between QC and RC regimes. Note that our interpretation of the numerical data is 
different from that in Ref. . 

To conclude, we have presented a theory of the critical properties of non-collinear quan- 
tum antiferromagnets in two dimensions. Our key assumption was on the validity of a 
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continuum description in SU{2) variables, which suppressed vortex excitations. However, 
we were then able to show that our results were consistent with earlier large and 
D = 2 + e studies. The quantum disordering transition was described by an anisotropic 
sigma-model for spin-1/2, bosonic spinon fields. All physical observables involve a collective 
mode of two spinous, and we computed explicit scaling forms for a variety of experimentally 
measurable quantities. Our results for Xu in the QC region are roughly consistent with 
recent numerical data on the spin-1/2 triangular antiferromagnet DTSf; this may be viewed 
as indirect evidence for the presence of deconfined spinous. However, numerical results also 
seem to indicate that the T range where QC behavior may be observed is rather narrow for 
this system. More detailed studies, especially in quantum-disordered non-collinear magnets 
will be quite useful. 
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